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ABSTRACT 

While emission and absorption lines of the 21 cm spin-flip transition of neutral Hy- 
drogen are intrinsically unpolarizcd, a magnetic field creates left- and right-handed 
polarized components through the Zeeman effect. Here we consider the resulting po- 
larization of the redshifted 21 cm background from the intergalactic medium before 
reionization. The polarization is detectable in regions with a strong gradient in the 
mean brightness temperature. In principle, this can open a new window on the evolu- 
tion of intergalactic magnetic fields. One possible approach is an extended integration 
of an individual target during this era, such as the Mpc-scale HH regions inferred to 
surround quasars at ^ ~ 6.5. The differential intensity between the two polarization 
states can be used as a probe of the magnetic field at the edge of the HII region. We 
estimate that the SKA could (ignoring systematics) detect B ^ 200 (10) /iG coherent 
over several kiloparsecs with an observational bandwidth of 100 (2) kHz. Beyond indi- 
vidual sources, the statistical properties of wide-field 21 cm polarization maps, such as 
the angular power spectrum, can be used to constrain the large-scale magnetic field. 
In this case, the SKA can detect B ^ 100/^G fields coherent over many megaparsecs. 
The magnetic field can be measured in any epoch over which the 21 cm background 
changes rapidly (for example because the ionized fraction or spin temperature change). 
Although the resulting constraints with SKA are relatively weak compared to theo- 
retical expectations, they nevertheless offer a unique direct probe of magnetic fields 
in the high-redshift universe. 

Key words: cosmology: theory — diffuse radiation — large scale structure — inter- 
galactic medium 



1 INTRODUCTION 

The next generation of low- frequency radio interferometers, 
such as LOFAR and the SKA, is expected to open up the 
possibility of observing neutral Hydrogen in the intergalactic 
medium (IGM) at redshifts ^ 10 when the universe transi- 
tions from neutral to fully ionized (e.g., Scott & Rees 1990; 
Madau et al. 1997; Zaldarriaga et al. 2004; Morales & Hewitt 
2004) through observations of the redshifted spin-flip tran- 
sition (with a rest wavelength of 21 cm) against the cosmic 
microwave background (CMB). This emission (during most 
of this era) and absorption (at high redshifts if the IGM is 
cold compared to CMB) is intrinsically unpolarized, so dis- 
cussions of the 21 cm background from the era of reionization 
have so far concentrated on measuring the total intensity or 
brightness temperature and its anisotropies. 

While the 21 cm emission is normally unpolarized, a 



magnetic field can change that. The Zeeman effect, in gen- 
eral, splits the signal into three parts: an unpolarized central 
component (tt component; Bolton & Wild 1957) bracketed in 
energy by two polarized components {a components). When 
viewed parallel to the magnetic field, the central component 
disappears and the two cr components are right- and left- 
handed circularly polarized, respectively. The difference in 
polarization states occurs because states split by the Zeeman 
effect have ±h momenta and, in order to conserve the total 
angular momentum, the two photons must also be polarized 
in opposite directions. Polarization of the 21 cm line has 
been successfully observed and used to determine the mag- 
netic field strength of the interstellar medium (Verschuur 
1968; Verschuur 1969) and that towards specific HI features 
(Verschuur 1995; Verschuur 1989) in the Milky Way, among 
other applications. 
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Observations of the polarized 21 cm background from 
high redshifts have two apphcations: (1) the direct deter- 
mination of the magnetic field strength in isolated regions 
containing neutral Hydrogen, such as the edges of HIT re- 
gions around quasars at z>6, and (2) the statistical deter- 
mination of large-scale magnetic field properties (such as its 
rms fluctuations) during the era of reionization. We shall see 
that the technique is most powerful when the mean bright- 
ness temperature of the 21cm line changes rapidly with fre- 
quency, for example because of sharp gradients in the ionized 
fraction or because of an overall heating of the neutral gas. 



This technique thus ofi^ers the potential to constrain 
the IGM magnetic field on a wide range of scales. At 
present, measurements of Faraday rotation against back- 
ground quasars constrains the IGM magnetic field to be 
< 10 nG at z ^ 0, if we assume a coherence length of ~ 1 
Mpc (see Kronberg 1994 for a review), but it could be larger 
if the field were tangled on small scales. Assuming flux con- 
servation, this would imply _B<([1 + z]/7)'^fiG, if the co- 
herence length remained 1 Mpc comoving. Within galaxies 
and galaxy clusters, field strengths are typically ~ 1-100/iG 
with coherence lengths of < 10 kpc (e.g., Clarke et al. 2001 
and references therein). The origin of these fields is uncer- 
tain. Early-universe processes can generate extremely weak 
seed fields on small scales (Quashnock et al. 1989). Fields 
can be amplified during the formation of galaxies through 
the dynamo process (Kulsrud 1999), but such a process will 
not work on the large scales of galaxy clusters (or the IGM) . 
Seed fields could also be generated during reionization itself 
(Gnedin, Ferrara, & Zweibel 2000). A more promising pos- 
sibility is for quasar or supernova winds to entrain the fields 
and carry them into the IGM (Kronberg ct al. 1999, 2001; 
Furlanetto & Locb 2001; Gopal-Krishna & Wiita 2001). 
Measurements of the field strength at 2 > 6 could offer in- 
valuable insight into these processes: because structure for- 
mation is much less advanced, galaxies and quasars should 
have polluted much less of the universe. However, placing 
such constraints promises to be difficult. For example, the 
number of bright sources available for Faraday rotation ob- 
servations is probably quite small, and the rotation may 
be dominated by the low-redshift universe because of the 
much larger path length. The polarized 21 cm background 
presents an alternate possibility; although we shall see that 
the constraints are fairly weak compared to most magnetic 
field generation scenarios, they are still the most promising 
direct measurements for the high-redshift universe. 



The discussion is organized as follows: in § 2, we discuss 
the polarized 21 cm signal and consider two separate appli- 
cations involving individual targets (§ 2.1) and large-scale 
statistics (§ 2.2). In the same sections, we also discuss the 
potential detectability of these signatures with upcoming 
low-frequency radio telescopes arrays such as the SKA. We 
conclude with a summary in § 3. Throughout the paper, we 
make use of the WMAP-favored LCDM cosmological model 
(Spergel et al. 2003). 



2 METHOD OF CALCULATION 

In the presence of a magnetic field, the Zeeman effect causes 
a level splitting of the 21 cm line with magnitude Sv = 
\i> — Vol (Lang 1999) such that 

eB 



= I'D ± 



47rmeC 



1420 MHz ± 14Hz 



B 



10 /iG 



(1) 



For typical magnetic fields of order ~ lOpG, the fractional 
separation of the two states is significantly smaller than the 
line width of the 21 cm emission (typically many kHz or 
more) . Thus, observing two separate peaks in the line profile 
is not possible. The polarization of these two states, how- 
ever, is a potentially measurable quantity if the magnetic 
field strength is such that the magnitude of the difference 
between left- and right-handed polarization is larger than 
the detector noise. 

We denote the magnetic field parallel to the line of sight 
as i3|| = n • B. We then write the total brightness temper- 
ature and the difference in the brightness temperature be- 
tween left- and right-handed polarization maps as: 



ATtot(i^) = ATLiu) + ATniiy) 
APdiff(i^) = \ATL{iy) - ATR(jy)\ 



(2) 



The difference term can be simplified because ATl{v) oc 
ATtot(t'o — Sv) while ATii(i^) oc ATtot(i'o + Sv), as the two 
states separated in frequency are also polarized in opposite 
directions. Thus 

^dATto 



APdiff(i^) 



= 2iAT,,,{v) 



n B 

10 tiG 



(3) 



where AT^^ati^) is the derivative of the emission profile with 
respect to the frequency (in Hz). Polarization observations 
have the advantage that, although 5v is extremely small, one 
can still constrain magnetic field with an observed band- 
width significantly larger than 5v, so long as the gradient 
AT^ot over that frequency range is large. 

During the era of reionization, there are two distinct 
possibilities for polarization observations: (a) individual tar- 
gets, such as HII regions surrounding z ^-^ Q quasars, and 
(b) maps of the high- 2; universe as a whole. We will first dis- 
cuss the extent to which magnetic fields in individual quasar 
regions can be measured with first-generation 21 cm inter- 
ferometers, such as LOFAR and the SKA and then discuss 
applications of statistical studies with wide-field maps of the 
21 cm background. 

2.1 Individual HII Regions 

It is now widely believed that the reionization process was 
both inhomogeneous and patchy on large scales (Barkana 
& Loeb 2001; Furlanetto, Zaldarriaga, & Hernquist 2004a). 
Luminous sources such as the first galaxies and quasars ion- 
ized their surroundings and these ionized patches grew and 
overlapped until they completely ionized the universe. These 
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ionized patches contribute to CMB temperature (Santos et 
al. 2003) and polarization anisotropies (Zaldarriaga 1997), 
cause fluctuations in the 21 cm background (Furlanetto, Zal- 
darriaga, & Hernquist 2004b), and correlate the two (Cooray 
2004a, b). The rcionization process was not instantaneous 
(Wyithe & Loeb 2003; Cen 2003; Chen et al. 2003; Haiman 
& Holder 2003) and could, in certain models, explain both 
the high electron scattering optical depth Tes = 0.17 ± 0.04 
measured through the largo angular scale bump in the 
CMB temperature-polarization cross-correlation (Kogut et 
al. 2003) and the rapidly evolving neutral fraction xh at 
z 6 towards Sloan Digital Sky Survey quasars (Fan et 
al. 2002). During the partially ionized phase, each high-z 
quasar should be surrounded by an HII region expanding 
into a largely neutral medium, with its size determined by 
the quasar luminosity and spectrum (Wyithe & Loeb 2004b, 
Mesinger & Haiman 2004). 

While 21 cm intensity observations of these sources 
could reveal the shape of these HH regions, especially in 
combination with Lya optical depth profiles (Wyithe & 
Loeb 2004b) , polaxization measurements can reveal the pres- 
ence of a magnetic field. To estimate the expected polariza- 
tion signal, wo only need the width of the edge of the HII 
region (which determines ATtotii') in eq. [3]). This is deter- 
mined by the mean free path A of an ionizing photon, which 
in turn depends on the photon energy (and hence the quasar 
spectrum) . Using a typical quasar spectrum, Wyithe & Loeb 
(2004b) estimate A ~ 1. 5x^^(1 + z)"'^ Mpc (proper). This 
corresponds to a frequency width of Ai^n ~ 2 kHz, much 
smaller than the resolution of any of the upcoming experi- 
ments. This is only a rough estimate given uncertainties in 
the exact shape of the quasar UV spectrum and the physical 
state of the IGM surrounding these quasars, but it suflices 
for simple estimates of the signal. 

Using the above case as an example, we estimate polar- 
ized signals of order APdiff ~ 0.3mK(_B|| /10/iG) if the ob- 
served bandwidth Ai^cb exceeds Aumi- Of course, only the 
transition region produces polarized flux, so the fractional 
polarization is ~ (At4iii/Ai^ch)(APdiff/ATtot) ~ 0.03% for 
l>^^^ = 0.1 MHz. Note that this assumes a field coherent 
across the several kiloparsec thickness of the transition re- 
gion; field reversals (either in frequency space or across the 
telescope beam) would cause the net polarization to exe- 
cute a random walk and weaken the limit. This weak signal 
presents a substantial challenge for an unambiguous detec- 
tion given the large number of systcmatics in these measure- 
ments. Nevertheless, with the rapid improvement in tech- 
nology and the unknown astrophysical environment (which 
could potentially increase the amplitude of polarized com- 
ponents), it is useful to estimate the level of magnetic fields 
detectable by future instruments. 

We let A^Ai^eh be the rms noise in each polarization com- 
ponent over a bandwidth Av^h- The minimum detectable 
magnetic field strength (parallel to the line of sight), at the 
signal-to-noise level of unity, is then 



Br (mg) 



V2N. 



(4) 



2.8AT4t(z/)[Az^ii/Az/ch] ' 
The factor \/2 accounts for the fact that the measurement 



is a difference between polarizations. The bandwidth factor 
appears because the signal contributes to only a fraction of 
the channel (we have assumed Ai/^h > Ai^n here). The root 
mean square noise contribution in a single SKA'^ channel is 



0.079mK 



^203MHz 



MHz y 



^eff/Ts, 



2 X 104m2K-l 
O.lMHz 30days\^/^ 



Avch 



tint 



(5) 



where AeB /Tsys captures the sensitivity of an array with a 
collecting area of ^efi and a system temperature Tsys, Ad 
is the beam size, and tint is the total integration time. The 
parameters are for the straw-man SKA, with revised values 
scaled from the estimate in Zaldarriaga et al. (2004). In one 
month of continuous integration, SKA can reach noise levels 
of ~ 80 /xK with these parameters, so one could potentially 
limit the magnetic field to a level of ~ 200 /xG at the one- 
sigma level. If we reduce the bandwidth to Avch = 10 (2) 
kHz centered on the edge of the HII region, the onc-sigma 
constraint is ^ 60 (9) /iG. While smaller bandwidths are 
desirable, resulting increase in systematics complicate such 
an observation. Narrow channels of order 2 kHz are also 
limited by the curvature of the HII region across the beam 
(which would cause the edge to move in frequency space). 
The latter could be avoided by careful binning of smaller 
pixels. 

Thus, for a reasonable bandwidth, the expected SKA 
limit is significantly higher than many of the theoretical ex- 
pectations at these redshifts. Models in which quasar out- 
fiows or starburst winds pollute the IGM usually have an 
average field strength of ~ 1//G in the magnetized regions 
(e.g. Furlanetto & Loeb 2001) and primordial magnetic fields 
are expected to be much weaker. Reaching a ~ 1 /iG de- 
tection level would require an order of magnitude increase 
increase in sensitivity relative to SKA and a better control 
of systematics such that observations can be made over a 
bandwidth of a few kHz level. However, it is important to 
note that strengths of ~ IfiG are only average levels. The 
fields could vary spatially within magnetized regions, espe- 
cially near galaxies and quasars. Because the edge of the HII 
regions are several kpc thick, this test would be sensitive to 
any strong fields tangled on these scales. 



2.2 Anisotropy Studies 

In addition to individual targets, one can also measure po- 
larization statistics with wide-field maps. This is analogous 
to proposed measurements of the 21 cm power spectrum 
(Zaldarriaga et al. 2004), and we take a similar approach to 
describing the polarized visibilities. 

In the case of polarization, we write the multipole mo- 
ments of a difference map of the two orthogonal states as 



diff _ / 
o-im — I 



dnAP^is(h)Yr''{A) 



(6) 



http://www.skatelescope.org for details 
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where APdifi(n) is now the difference in maps when obser- 
vations have integrated over a finite bandwidth such that 

APdiff (n) = 2.8 j drAT^,,{nr,r) [h ■ B(nr,r)] W^r) . (7) 

Here the magnetic field is measured in units of /iG, r is 
the comoving distance, and W^{r) describes the experimen- 
tal response, assumed to peak at a distance ro. The window 
function thus picks out a particular redshift interval because 
the 21 cm feature is a spectral line. We take the response 
function to be a Gaussian of width Avch- The frequency 
derivative can also be transformed to a function of redshift 
via dAT{z) /du = —vo/v^dAT(z)/dz, where u is the ob- 
served frequency, v — U{){1 + z)^^ . 

The maximum magnitude is set by the overall strength 
of the 21 cm intensity. Following Zaldarriaga et al. (2004), 
the mean brightness temperature of the 21 cm line is 



23.2 mk 



0.02 



xh{z) 



Ts — Tch 



Ts 



(8) 



where Tcmb = 2.73(l + .z) K is the CMB temperature at red- 
shift 2 and Ts is the gas spin temperature. Henceforth for 
simplicity we will neglect spatial fluctuations in Ts and xh 
(though we will allow both to be functions of redshift). The 
angular power spectrum of the polarization difference map 
includes fluctuations both in the neutral Hydrogen density 
field (which are captured by eq. and spatial fluctuations 
of the magnetic field. With our approximation, fiuctuations 
in the density cause 5T{n) = AT^'^'"^{n)[l + Sg{n)], where 
uh = XHng{l + Sg), the mean baryonic density of the uni- 
verse is fig, and the mean neutral fraction is xh- They can 
thus be encapsulated in Pgg{k), the spatial power spectrum 
of gas density fiuctuations (with respect to the mean den- 
sity). The remaining term depends on PBs(fc), the spatial 
power spectrum of the magnetic field (which we will mea- 
sure in absolute units of nG^ Mpc"''). The relation between 
the two power spectra is currently unknown. It is possible 
that the coherence scale of the magnetic field is larger than 
that of the density field (this is similar to the the velocity 
field, which has a larger correlation length than the density 
field). On the other hand, the opposite situation could apply 
if magnetic fields are highly tangled. 

In either case, the power spectrum may then be written 



^diff / dr 



2.8- 



dAr 



1^0 



dz 



/(J^;r) 1^.(0' (9) 



in the Limber approximation appropriate for a fiat sky (Lim- 
ber 1954). Here the mode coupling integral / is a convolution 
between fiuctuations in the gas density field and the mag- 
netic field, weighted by a geometric term involving the line 
of sight angle. 



/(fc) = 



djj. 



PgAkV2)PBB{kyi), (10) 




Redshift 

Figure 1. Mean brightness temperature as a function of the red- 
shift for two different models of the reionization history. In model 
A (dashed line), we have ^ = 20, /* = 0.05, and /heat = 1- Here 
heating occurs early and ^7^2lcm ^i-g^ces xh for z < 15. Model 
B (solid line) has = 6, /* = 0.03, and /heat = 0.01. Heating 
is much less rapid in this case, but A2"2lcm g^jjj ij-q^^^^ xjj for 
z<ll. 



the mean brightness temperature and fiuctuations in the 
line-of-sight magnetic field. The same formula can also be 
derived as an overall modulation of the mean brightness 
temperature by fiuctuations in the magnetic field and, in 
this sense, the derivation is analogous to calculations in the 
literature (e.g., Cooray 2001) of the Ostriker & Vishniac 
(1987) effect, which is a modulation of the density field by 
the velocity fiuctuations. 

In the limit in which small-scale fluctuations in the 
brightness temperature (from density structure) modu- 
late a more coherent magnetic field, ^2 = |k — ki| ~ 
k such that j/i ^ 1 and y2 1. Thus, I{k) — > 

Pgg{k) J kldki/2TT^PBB{ki) J^^ dfi{l ~ fi^)/4. The integral 
over jj, simplifies to 1/3 such that I{k) 
where 



Pa9ik)BL.ir)/3, 



By 



k^dk 
27r2 



PBB[k). 



(11) 



Thus, if the magnetic field is spatially uniform on larger 
scales than the density fiuctuations, we have 



^diff _ / dr 



dv 



(12) 



In this case, the factor l/3-B'^)rms can also be understood 
as a random averaging when the two line of sight angles are 
nearly parallel (i.e. the small angle approximation): 

1 



[n.B(n)][n'-B(r2)]) 



(13) 



where j/i — ki/k and y2 ~ k2/k = -y^l — 2/ij/i + yf. This 
form for the power spectrum occurs because equation (7) 
presents a convolution between the frequency derivative of 



We will now restrict ourselves to the simple limiting 
case of eq. 11211 appropriate for large-scale magnetic fields. 
To proceed, we need a model for dAT/dz. This has three 
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components: the cosmological density evolution, xh, and Ts- 
We compute xh following Furlanetto et al. (2004a), who as- 
sume that the ionized fraction is proportional to the collapse 
fraction (with a constant C^, part of which is determined by 
the star formation efficiency /*). The spin temperature evo- 
lution is more subtle. At these redshifts, Ts is coupled to the 
gas kinetic temperature by Lya photons (Field 1959; Madau 
et al. 1997). We follow this background in the same way as 
Ciardi & Madau (2003) , assuming that five photons are pro- 
duced near Lya for every ionizing photon, as appropriate 
for Population II stars with a constant star formation (see 
their Fig. 1). We assume that the IGM cools adiabatically 
after decoupling from the CMB. We then heat it through 
X-rays from supernovae inside the ionizing sources. We as- 
sume one supernova occurs per 100 Mq of stars and that 
each has 10^^ ergs. We further assume that a fraction /heat 
of the supernova energy escapes into X-rays and uniformly 
heats the IGM. Fig. 1 illustrates two models. Model A is 
a 'maximal heating' scenario in which /\y2icm j-jggg j-^pidly 
and decreases later only through reionization. Model B is a 
'minimal heating' scenario in which the IGM remains cool 
for a long period but still emits well before reionization oc- 
curs. Neither of these models can be taken as more than a 
representative possibility, as they neglect a number of im- 
portant effects (such as shock heating in the IGM, other 
X-ray sources, and a rigorous treatment of the Lya back- 
ground), but they suffice for the purposes of estimating the 
brightness temperature gradients for our calculation. 

In Fig. 2, we show the polarization and temperature 
angular power spectra for model B at z = 9. In these cal- 
culations, we describe the density power spectrum with the 
linear density field, but scaled by a bias factor appropriate 
to the distribution of neutral Hydrogen (Mo et al. 1997). 
Specifically, we follow the methods of Santos et al. (2003, 
2004), who use an approach similar to halo models of large 
scale structure (Cooray & Sheth 2002). We also show noise 
curves for the SKA and for an 'extended SKA'. The po- 
larized signal is proportional to the brightness temperature 
gradient. We have chosen a point at which this is relatively 
shallow in Fig. 1; if the rapid temperature evolution of model 
B at 2 ~ 15-25 is accessible to observations, the polarized 
signal could be a few times larger. The signal for the late 
stages of reionization in model A would be comparable to 
that shown in Fig. 2. Moreover, we have assumed a uni- 
form Xh; spatial variations in the ionized fraction (or in Ts) 
could also increase the signal by a factor of a few. On the 
other hand, we have neglected possible correlations between 
the density and the magnetic field, and we have neglected 
small-scale structure in B. Our estimate can therefore only 
be taken as a guide to the expected order of magnitude of 
the signal. 

We can estimate the smallest detectable Brms using the 
Fisher matrix approach. Since only a single parameter is 
involved, we can write the limit on the magnetic field as 
simply 



where 



I 



1 (dCf^ 



(14) 



10- 
10" 
10" 
10" 
^.10- 

E 

^ 10" 
10" 
10" 

io" 





- 







- 






lOOnG 












10[iG 




Polarization 




■ ^ 







10" 



10' 



10^ 

/ 



10^ 



10' 



Figure 2. The rms brightness temperature fluctuations of the 
21 cm background at 2 = 9 (i^obs = 142 MHz), where ATrms = 
+ l)C;/27r. The top solid curve (labeled 'Temp.') shows the 
total temperature fluctuations, while the long-dashed curves show 
the fluctuations in the polarization difi'erence assuming a large 
scale magnetic field with -Brms = 10 and 100/iG. Both assume 
model B of Fig. 1. The 'SKA' curve shows the instrumental noise, 
unbinned in the multipole space, for a one month integration with 
Ai/^h = 0.5MHz. The bottom (linear) curve assumes an order of 
magnitude improvement in sensitivity over the SKA, a year-long 
integration, and a total bandwidth of 8 MHz. 



{21 + l)/eky 



Cnoisc \ 



(15) 



with /sky the fraction of sky observed and (7j"°'=<= the instru- 
mental noise. We make the null hypothesis Cf'^ = and 
estimate the one-sigma constraint on Brms for the two dif- 
ferent noise curves shown in Fig. 2. For the standard SKA 
noise parameters with a one month observation), we find 
Brms > 224(/sky/0.1)"°-^ IJ.G IS detectable (at la). With 
an order of magnitude improved sensitivity and a year-long 
observation, the limit improves to Brms > 7.5(/sky/0.1)~'''^ 
nG. As with targeted observations of HII regions, it is thus 
difficult to probe the ^G fields expected from many models. 
Long integrations with instruments beyond the SKA will be 
required to reach the necessary limits. 



3 SUMMARY 

While emission (and absorption) from the 21 cm spin-flip 
transition of neutral Hydrogen is intrinsically unpolarized, 
a magnetic field creates left- and right-handed polarization 
through the Zeeman effect. For individual targets of neutral 
Hydrogen during the era of reionization, such as Mpc-scale 
HII regions surrounding quasars at ^ 6, the difference in 
intensity between these two polarization states can be used 
as a probe of the magnetic field on the scale of the edge of the 
HII region. For example, observations of the known SDSS 
quasars at z ~ 6.5 with the SKA can constrain magnetic 
fields to a level of 200 juG (at the la level) in an integration 
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of one month over a bandwidth of 100 kHz. If a few kHz 
bandwidth observations can be made with systematics well 
under control, the limit improves to ~ 10 fiG, again at the 
Icr level. The expected fields have B < 3/iG so one could 
potentially reach an interesting constraint through a long 
integration or an even larger instrument. In addition, these 
measurements are sensitive to fields with scales of several 
kpc, which could be much larger than the mean value. 

Beyond individual sources, the statistical properties of 
21 cm polarization in wide-field maps, such as the angular 
power spectrum of fluctuations, can also be used to constrain 
large scale fluctuations of the magnetic fleld. This requires 
the mean brightness temperature of the 21 cm background 
to vary rapidly as a function of the frequency. Thus polariza- 
tion detections are limited to epochs over which the underly- 
ing properties of the 21 cm background, such as the ionized 
fraction and spin temperature, change rapidly. Using simple 
models for the evolution of the mean brightness tempera- 
ture, we have estimated the statistical constraints that can 
be placed on a magnetic field coherent over scales larger than 
that of the density field for such a measurement. The SKA 
can detect a large-scale field with Brms > 224(/aky/0.1) "'^ 
jiG. A future instrument with an order of magnitude better 
sensitivity can improve the limit to Brms > 7.5(/sky/0.1)~°'^ 

with a year-long observation. Constraints of this order 
would be desirable given that galaxy clusters at 2; < 1 are 
known to have firms ~ 10-100/iG (e.g., Clarke et al. 2001). 

Finally, we emphasize that, although the limits obtain- 
able from the 21 cm background are weak compared to the- 
oretical expectations, obtaining other direct constraints on 
the magnetic fleld during the era of reionization promises 
to be equally challenging. This technique may offer the first 
direct limits on high-redshift magnetic fields. 
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